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Abstract 



' We determine exactly the short-distance leading behavior of the density correlation functions of 

^ a two-dimensional two-component charge-symmetric Coulomb gas composed of point particles, in 

the whole regime of stability where the coulombic coupling (3 <2. More generally, we compute the 

' short-distance behavior of the effective interaction potential between two external arbitrary charges 
^ . 

■ Qi and Q2 immersed in the plasma, for /3|(5i| < 2 and /3|(52| < 2. We also find the short-distance 
g. 

I , asymptotics of the density profiles near a single external charge Q immersed in the plasma for 

g; /3|Q|<2. 
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I. INTRODUCTION 



The two-dimensional charge-symmetric two-component plasma is a two-dimensional 
Coulomb fluid composed of two species of point particles with charges ±1. Each pair of 
particles (charges q and q') at a distance r from each other interact via the logarithmic 
Coulomb potential in two dimensions 

T 

v{r) = —qq'\n— (1.1) 

where L is an arbitrary length scale. We describe the system with classical (i.e. non quantum) 
statistical mechanics at a reduced inverse temperature (coulombic coupling) j3. Provided 
f3 < 2 the system of point particles is stable against the collapse of particles of opposite 
sign. We will consider here only the regime (3 < 2. To study the regime beyond /3 = 2 it is 
necessary to consider an additional repulsive short-range interaction to avoid the collapse. 

The two-component plasma is mappable into the quantum sine-Gordon model Q], which 
in, two dimensions, is an integrable model For (3 < 2, using results from the 

sine-Gordon model, the bulk thermodynamic quantities of the two-component plasma have 
been computed exactly Isf. Also, using recent results for the form factors of the sine-Gordon 

un 

model p, Iq], the large-distance behavior of the particle correlation functions have been 
determined exactly In Ref. Q|, the large-distance behavior of the mean- force effective 
potential between two external charges immersed in the plasma was determined. Also in 
that same work ^1^, the large-distance behavior of the density proflles around an external 
charge in the plasma was computed (see also 111'). 

n n n 

The work presented here is complementary to the ones of Refs. H M UJ. We are 
interested here in the exact short-distance expansion of the effective interaction energy 
between two external charges Qi and Q2 immersed in the plasma. We will extensively 
study the cases Q2 = ^1 and Qi arbitrary and determine the short-distance expansion of 
the density proflle of around a single external charge in the plasma. Also we will study the 
case Qi = 1 and Q2 = ^1 which allows us to determine the short-distance expansion of the 
density correlation functions. Our main tool will be the operator product expansion of the 
exponential flelds of the sine-Gordon model. 

This paper is organized as follows. In Sec. |n] we will recall some basics facts about 
the exact solution of the two-component plasma and the sine-Gordon model needed for 
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our purposes. In Sec. IIIII follows a general discussion on the short-distance expansion of 
the mean-force potential interaction between two arbitrary external charges Qi and Q2 
immersed in the plasma, which is obtained by means of the operator product expansion. 
In Sec. IIVI we will consider the case of charge-like particles Q1Q2 > and in Sec. |V] the 
case of particles with opposite signs Q1Q2 < 0. Both sections IIVI and |3 have a similar 
organization, although some technical details are different. First, we specialize in the case 
where = IQ2I (identical or opposite external charges) and we study the mean- force 
potential between these two external charges. Then, we consider the case where (^2 = =tl is 
a charge of the plasma and Qi = Q > is arbitrary and we study the density profile around 
the charge Q. Then, we consider the case Q2 = ±1 and Qi = 1 and determine the short- 
distance behavior of the density correlation functions of the plasma. We test our general 
results valid for arbitrary jS < 2 against the known results for the correlation functions when 
(3 = 2. In the appendices are presented some expansions for particular values of /?, Qi 
a,.d of the expectafo,. vahjes of exponential fields and of so,ne Conlon>b-type .ntegrals 



(Dotsenko-Fateev integrals 



12, ll^) that are needed in the main text. 



II. THE TWO-COMPONENT PLASMA AND THE SINE-GORDON MODEL 

The statistical mechanics of the two-dimensional two-component plasma are best worked 
out in the grand canonical ensemble with fugacity C (with dimensions length'"^). Writing 
down the Boltzmann factor of the system one can notice that only the combination z = C,L^^'^ 
is relevant. We shall still call z the fugacity, although it should be noticed that its dimensions 

are length^ ^"^'"^ . 

Doing a Hubbard- Stratonovich transformation, the grand canonical partition function S 
of the plasma can be written as the partition function of the quantum sine-Gordon model 



„ _ / r'0exp[-5(z)] 
/ P(/>exp[-^(0)] 



with 



and 



S{z) = - 



IQ-K 



+ 2z cos(60) 



13/ A. 



(2.1) 



(2.2) 



(2.3) 
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Under this mapping, the bulk density and two-body densities of charges q = ±1 and q' = ±1 



are expressible as 



nl = z(e^^^^) (2.4) 



and 

nf^,{\r - r'l) = (^e''"^<^>ir)^^bg'Hr')^ ^2.5) 

where the averages are taken with respect to the sine-Gordon action (j2.2|) . 

In the sine-Gordon mapping the fugacity z is renormalized by the (diverging) term 
g/3f(o)/2 Thus, to make a complete connection between the sine-Gordon model and 

the statistical mechanics problem of the two-component plasma it is necessary to specify 
the proper renormalization scheme for the fugacity. In Ref. jy|, Zamolodchikov considered 
the sine-Gordon model as the Gaussian conformal field (with action S{0)) perturbed by 
the cos{b(f)) operator. Under this scheme, he used the conformal normalization, where the 
exponential fields e*^*^ are normalized such that 

(^(,ib<i,i^i) . . . g»b0{x„)g-i6<A(yi) . . . e-^'"^(y"))^^o = ni<jo<Jx» - Xjl^ly, - Yjl^ ^^^^^ 



which is clearly the appropriate normalization for the Coulomb gas system in view of rela- 
tions (Q and (Q. 

The specific dimensionless grand potential of the sine-Gordon model was computed in 
Ref. yi in terms of the mass of the soliton of the sine-Gordon model. On the other hand 
the relationship between the fugacity and the mass of the soliton under the conformal nor- 
malization was found in Ref. With these results the density-fugacity relationship and 
full thermodynamics of the Coulomb plasma can be obtained 9] . From the fugacity-density 
relationship from Ref. an explicit formula for the expectation value of the exponential 
field, which we will need in the following, can be obtained (see also j^) 

^ ^ Kim)) W-f)) (4 - /5)7(/5/4) ^'-'^ 

where we use the notation 7(x) = r(x)/r(l — x), with T{x) the Euler gamma function, and 
e = /?/(4-/?). 
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III. THE INTERACTION POTENTIAL BETWEEN TWO EXTERNAL 
CHARGES AND THE OPERATOR PRODUCT EXPANSION 



Let us consider that two external point charges Qi and Q2 are immersed in the plasma 
and located at the origin and at r respectively. To avoid the collapse of the external charges 
with the charges of the plasma of opposite sign, we suppose that l3\Qi\ < 2 and /3\Q2\ < 2. 
We define the effective interaction mean-force potential EQ^Q^{r) between the charges Qi 
and Q2 as in [l^ by 

S[Qi,0;Q2,r]/S 



exp[-/3EQ,Q2(r)] 



(3.1) 



(H[gi]/S)(S[Q2]/H) 

where S[Qi,0;Q2;r] is the grand partition function of the plasma in the presence of the 
charge Qi at the origin and of the charge Q2 at r (with the direct interaction —QiQ2\iir 
between the external charges included in the Hamiltonian) and 'E\Q] is the grand partition 
function of the plasma in the presence of a single external charge Q. As shown in Ref. lol |. 
the effective interaction can be expressed in terms of expectation values of exponentials fields 
of the sine-Gordon model 

|/gi6Qi</.(0)gi6Q2</>(r)^ 



exp[-/5^QiQ2(r)] 



(3.2) 



The denominator in the preceding expression can be computed since an exact expression 
for the expectation value _(e^''^ was conjectured in Ref. Q], a conjecture later supported by 



several other studies 



15 



16 



id, 



with 



hiQ) 



~ dt 
T 



smh.\2QbH) 



2,2 ~2t 



(3.3) 



(3.4) 



_2sinh(62t) sinh(t) cosh[(l - 1)^)1] 
This expression is valid for [3\Q\ < 2, otherwise the integral hiQ) diverges. Notice that this 
is the same stability condition necessary to avoid the collapse of an external charge Q with 
the counterions of the plasma (ions with charge of opposite sign than Q). 

On the other hand, the short-distance expansion of (^e^^Q-i't'W ^'^''Q^'I'i^)'^ can be obtained by 
means of the operator product expansion [l8^ 



gibQi</.(O)gjbQ20(r) 



which reads 



19| 



n=+oo 



ib{Qi+Q2+n) 



+ 



(3.5) 



where the dots denote subdominant contributions from descendant fields (hke 
(9(/>)^(9(/))^ e*^*^^^'^^'^'^'''^ and so on), which will not be needed for our purposes. The co- 
efficients of the operator product expansion have the following form j3] 

(jnfl |nl^/3Qi(32+n/3(Qi+Q2)+2ln|(l-|)+n2/3/2p,0 (^'ijA-P\ g^ 

Q1Q2 ^ ' QiQ2^ ' V ■ / 

where Jq^q^ admit a power series expansion of the form 

00 

/S.(^) = E/^'°(^l'^2)x^ (3.7) 

k=0 

The coefficients of the first leading terms of the expansion ()3.7|) are 

C{Qi,Q2) = 1 (3.8) 

/o"'°(Qi, Q2) = UPQi, PQ2, P) forn > (3.9) 

fo'\QuQ2) = J\n\{-f3Qi,-PQ2,f3) forn<0 (3.10) 

where 

Jn(/9Qi,/5Q2,/3) = - / l[d^Xkl[\xkr^\l-Xkf'^' II \xk-xj\^. (3.11) 

■ k=l k=l l<k<j<n 

These Coulomb-type integrals have been computed exactly by Dotsenko and Fateev lisj 

n 

jniPQi, (3Q2, = vr"7(/3/4)-"n7(W4) x (3.12) 

k=l 

II7 (^1 + ^ + fc^J 7 (^1 + ^ + fc^J 7 (^-1 - 2 + Q2)-{n-l + k)-j . 
The coefficient for n = of the first sub leading term {k = 1) in (|3.7|) reads 

/r(Qi,Q2) = J(/3Qi,/3Q2,/3) (3.13) 

where 

J{f3Q,,(3Q2,f3) = / d'xd'y ^ , ' ' ' (3.14) 



which can be expressed in terms of hypergeometric functions 3F2 at unity j20|,|21[. In Sec. IVl 
and in the appendix |Bl we will compute a special case of this integral when Qi = —Q2- 
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Putting all these expressions together, the first terms of the short- distance expansion of 
the Boltzmann factor of EQ^Q^{r) are 

lpih(Qi+Q2)'i>\ 

exp[-/J£«.<,.(r)l = r"** (3.15) 

n=l ^ ' ^ ' 

n=— oo ^ ' ^ ' 

/pjb(Qi+Q2)0\ 

+ (e^Q..)(e-J^) -^-'"^"^"^-'^ + + 0(r^«^«-^) . 

The 0(r^*5^'^2"'"'*) ^g^jxi in the last line comes from a descendant field contribution which is 
explicitly computed in Ref. Q], but that will not be needed in the following. 

If (3Q2 and (3 are small enough the dominant contribution is given by the term n = 
in (jS3I) 

eM-pE^.Q.{r)]r.^ ^^J,^ r^^^^^ (3.16) 

which is an expected result: the two charges interact at short distances with the bare 
Coulomb potential — QiQ2lnr. Notice however that depending on the values of [3Qi, PQ2 
or f3, some of the following terms in the expansion ()3.15|) can become dominant over the 
n = term. We will explore this possibility in the following sections. 



IV. THE CHARGE-LIKE CASE: Qi > AND Q2 > 

Let us consider that Qi and Q2 have the same sign. In the general small-r asymp- 
totics ()3.15j) . the n = —1 term will be dominant over the n = term when (3{Qi + Q2) > 2, 
a condition which will take place if Qi + Q2 is positive and sufficiently large. On the other 
hand the term n = 1 will be dominant over the n = term if f3{Qi + Q2) < —2. As we see, 
depending on the sign of Qi + Q2, the terms of the series for n < (if Qi + Q2 > 0) oi n > 
(if Qi + Q2 < 0) will become dominant. 

Let us suppose, without loss of generality, that Qi > and Q2 > 0. The term —n {n > 0) 
will be the dominant term of the series ()3.15p if 2 + (ra — < (3{Qi + Q2) < 2 + n[3. Thus, 
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the dominant term of the small-r expansion of exp[— /3i?Q^Q2(r)] is 



2 + - < /3(Qi + Qa) < 2 + n = 1, 2, . . . 



This relation is vahd provided we are still in the region of stability of the system, thus 
PQi < 2, PQ2 < 2 and /? < 2. 

There is an interesting relation between the changes of the dominant small-r behavior of 
the effective potential shown in ()4.ip and the corresponding coefficient expressed in terms 
of the Dotsenko-Fateev integrals ()3.1ip . jn{—PQi,—PQ2,P) expressed as an integral in 
equation can be interpreted as the partition function of a Coulomb gas with the 

external charges Qi at and Q2 at 1 and n charges —1. Thus, the behavior of EQ^Q^{r) 
when 2 + — < l3{Qi + Q2) < 2 + n(3 is dominated by the approach of n charges —1 
to the positive external charges Qi and Q2- This confirms and generalizes the analysis of 



Hansen and Viot 



221. 



The integral (jTTTll defining jni-PQi, -l3Q2, P) converges at — >■ if f3Qi < 2, at 
Xjfc ^ 1 if PQ2 < 2 and at \xk\ — > 00 if j3{Qi + Q2) > 2 + (n — On the other hand, 
|/g«b(Qi+Q2-n)0^^ converges if P{Qi + Q2) < 2 + njS. These are precisely the same conditions 
for the corresponding term in the small-r expansion ()3.15p to become dominant. Thus 
when the term —n of the general operator product expansion ()3.15|) becomes dominant, its 
corresponding coefficient in ()4.H) is properly defined. 

When the condition + Q2) > 2 + (n — 1)[3 is not satisfied, the integral defining 

in in equation (jH.llj) is not convergent, but the corresponding term in ()3.15|1 is not the 
dominant one. To evaluate this subdominant contribution one has to resort to the analytic 
continuation of j„ provided by the expression ()3.12j) in terms of Gamma functions. 

At j3{Qi + Q2) =2 + {n — 1)13, the term —n and the term —n + 1 of ()3.15p are of the same 
order, in 7-WiQ2-n{n-i)/3/2^ they should be added to obtain the correct leading behavior. 
The coefficient of the term —n exhibits a pole of first order due to jn{—l3Qi, —l3Q2,l3), but 
the coefficient of the term —n + 1 of the expansion ()3.15j) also exhibits a pole of first order 
due to (e^^(«i+Q2-"+i)'^). Indeed, if we write p{Qi + Q2) = 2 + {n - + pe, we have Q 

Uib{Q,+Q2-n+l)^ /e»b{^-i)\^ (4 2) 
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and 



Jn(-/9Qi,-/3Q2,/5) = vrj„,-i(-/3Qi,-/5Q2,/5) 



7(f)7(f + f) 
7(f + f) 



7(1 - /?gi + ^ ^^ )7Wi + ^ ^ - y )7( y) 
~ Jn-i(-/3Qi,-/?Q2,/3)|^. 



(4.3) 



Both coefficients turn out to have the same residue in absolute value but with opposite 
signs, and thus cancel each other, giving a finite value, and an additional 7'/3QiQ2-n(n-i)/3/2 ^^j. 
contribution to the small-r expansion. This will be illustrated in several particular cases in 
the following sections. 

A. Effective mean-force potential between two identical external charges 

In this section we consider two identical external charges immersed in the plasma Qi = 
Q2 = Q > 0. Using equation (j4.ip . the first few changes in the small-r asymptotics of the 
effective potential read 

lAb{2Q-l)(t>\ 



-PEqq{t) 



^gi6(2Q-2)0^ 



< /3g < 1 

l< (3Q <l+ (3/2 



1 + 



(n-l)/3 



l + f</?g<l + /3;/?<l 



</?g<l + f ;/3<^;n = 3,4,... 

(4.4) 



Notice that in order to satisfy the relation I3Q > 1 + {n — l)(3/2 and the stability condition 
PQ < 2, the coulombic coupling must be small enough: /? < 2/(n — 1). Thus, as P becomes 
larger there are few changes in the small-r behavior of the effective potential. 

B. Colon density profile around an single external charge Q > 



Let us consider now a single external charge Q > at the origin immersed in the plasma. 
The density of the particles with charge -|-1 (colons) of the plasma is given by |lOl 111] 



n^(r) = z 



l^^ibQ<i>{SS) ^ib<i>{r)\^ 



(4.5) 



or equivalently, using relation ()2.4p . 

"+(^) = < \e-^Q^)(e^.^) (4-6) 

Thus, using the expansion ()3.15p with Qi = Q and Q2 = +1 we can deduce the small-r 
asymptotics of n+(r). 

To ensure the stabihty of the system we require that f3Q < 2. This condition reduces the 
number of changes in the small-r behavior (j4.ip to two. Writing explicitly the ji coefficient 
we find that the dominant behavior of the colon density is 

n+(r) ~ <^ ^ (4.7) 



7(l-f)7(l-f)7(-l + ^)r^-^ 2-/?</?Q<2 



Where (^e^^(Q+'^)'l''^ and (e*'"^'^) can be computed from equation (j3.3|) . This confirms and 
completes the analysis of Ref. where the change of behavior at l3Q = 2 — /3 of the colon 
density was identified and interpreted as a "precursor" of the counterion condensation, in 
the sense that for jSQ > 2 — (3, the effective potential (times (3) of a colon and the charge 
Q behaves at short distances as the bare Coulomb potential with a reduced charge: — (2 — 
/3) Inr, instead of —jSQlnr. The charge is reduced because of the counterion condensation. 
With the present work, we have also an explicit expression for the coefficient multiplying 

It is also interesting to study the transition at f3Q = 2 — (3 and at (3Q = 2. At j3Q = 2-/3, 
both terms r^'^ and r^~^ shown in equation ()4.7|) become of the same order and should be 
added to obtain the correct dominant behavior. 

Using the expansions of the expectation values of exponentials field ()A2jl and ()A6j) com- 
puted in the appendix |X1 we have 



■X 



2-^(0 + 1) 



Ri{b ) + hiT^ - 1) - 1 n In 



^262 ^ 4-/5 V7(^^) 



+ 0([/3(Q + 1)-2]2)| 



(4.8) 



where I^Q) = dIb{Q)dQ and Riip'^) are given in equations ()A3|) and ()A7|) of the ap- 
pendix IXj). On the other hand the 7( — 1 + [3{Q + l)/2) from ji also exhibit a pole at 
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j3Q = 2 — P, but with opposite residue: 



^2:7 1 1 - ^ 1 7 (1 - 2 1 ^ 



27rz 



2-(3{Q + l] 



-1 



[2C + ^(1 - /5/2) + ^iP/2)] + OmQ + 1) - 2]') 

(4.9) 

'(/'(I) is the Euler constant. 



where ^p{x) = T'(x)/T{x) is the digamma function and C = 
Thus adding both contribution gives a finite resuh. Using ()4.8|1 . ()4.9j) and expanding r^^ 
r^-f^il + [(3{Q + 1) - 2] Inr + 0{[f3{Q + 1) - 2]^) finally yields, for (3Q = 2-f3, 



In 



4-/3 



1 



262 



^-2 



L - /?/2) + V^(/3/2)] 

(4.10) 

Notice that the power of z in the logarithm, 

^2/(4-/3)^ is the appropriate one to have a 

dimensionless argument in the logarithm. 

When PQ = 2, the contribution from the term n = —2 of in f-PQ-P+'^^ jg of the 

same order as the one for n = —1, in r^"'^. The coefficients of both terms exhibit a pole 
with opposite residues when (3Q = 2. A similar mechanism as before applies, adding both 
contributions gives a finite result. The term n = —2 of the expansion of n_^{r) reads 



(pibiQ-l) 



with 



z r 



-j2(-/3g,-/3,/3) 



(4.11) 



PQ 



PQ , P, 



PQ , P. 



j2(-/5Q,-/3,/5)=7r^7(l-T) 7(1-^)7(-1 + ^)7(-1 + ^ + t)7(1-^ + t) (4-12) 



4- 



4' 



When PQ ^ 2, we have 



-j2i-PQ,-P,P) 



2ti 



PQ^2 z{PQ - 2) 



1 + 



PQ-2 

P 



Ri{b')+n{—-l)-2^\n 



TCZ 



•262 ' ^ 7(62 
while from the term n = — 1 we have a contribution from the ji coefficient: 

PiQ + l)\ 2 



+ 0i{PQ-2f)\ (4.13) 



■1 + 



/3Q^2 2- PQ 



-2C-ij{l-P/2)-^{P/2). 

(4.14) 
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Expanding r-^'?-^+'^ = r^-z^fl + (2 - pQ)\nr + 0((2 - PQ)'^)], and adding all contribu- 
tions ()4.13p and ()4.14|) . we actually find that at /3Q = 2, the small-r behavior of n+(r) is 
exactly the same as the one for j3Q = 2 — (3 given by equation (|4.1U|) . 



(2) 

C. Density correlation function n_^^{r 



The two-body density correlation function between particles of the same sign (say posi- 
tive), can be expressed as 

^gi60(O)gib(^(r)<^ 



b \2 



(4.15) 



Thus, using the results of the last section with Q = 1 allow us to find the small-r behavior 
of the correlation function: 



2TTz'^n^, r <{ — In 



7(1/4) 



2/3 



+ i?i(l/4) + /' 2(1) + 21n2 



p < 1 
P = l 



7iz^n\-^{l - If-fiP - 1) r^~t^ = nz^{e""^)'y{l - lf-f{p - 1) r^-f^ 1< /3 < 2 

(4.16) 

A numerical estimate of the constant intervening in the case P = 1 is Ri{l/4) + 7^^2(1) — 
-1.73246. 

Our result ()4.16|) confirms the change of behavior of the correlation function (r) at 



P = 1, predicted by Hansen and Viot uM- Furthermore, we obtained an explicit analytical 
expression for the coefficients multipliying (for P < 1) and r^~^ (for P > 1). 

The limit P ^ 2 cannot be obtained from the results of the previous section because 
there is an additional divergence due to n^_^. When /? — * 2, the term n = — 1 of the general 
series ()3.15p . in r^~^ is of the same order as the term n = —2. This term (n = —2) gives a 
contribution to n^^(r) equal to 



vrW-^(l-^)^(l-|)7(§-l)7(-l + f ) 



4^2^4^4-2/3 p 3^ 

-7(l-93^(-l + ^). (4.17) 



(2-/3)2 4' " ■ 4 

Notice that the coefficients of r^"^ in ()4.16p and of r'^"^/? f|4.l7j) exhibit now a pole of 
second order when /? — 2, as opposed to the situation of the previous section, when the pole 
was of order one. Indeed, 

-47(l-f)M-l + ¥) 4 



{P-2f 



/3-2 {p - 2y 



+ 0{P-2) 



(4.18) 
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while, using the expansion ()A9|1 of {e '''), and expanding the gamma functions, 



(e-)7(l-5)^7(/3-l)=^,^_,). 



l + (^-2)[C + ln(7rz)] 



(3 - 2V 

+ ^ [C + H7rz)][l + C + \n{nz)] 



+ 0iP-2) (4.19) 



Notice that now the poles do not cancel each other. Expanding ^ = 1 + [2 — f3)\nr + 
M)!(inr)2 + 0((/?-2)3) and r^-^f^ = 1 + 2(2-/?) Inr + 2(2 - /?)2(lnr)2 + 0((/? - 2)^), and 
putting together all contributions from ()4.18|) and ()4.19p . we finally find 



p-2 



[C + ln(7r2)] 



- [ln(7rzr) + C]^ + 2[C + \n{nz)]^ + C + \n{nz) + 0{/3 - 2) 



(4.20) 



In the preceding expression we recognize the expansion ()A10|) of {e^^'^)'^ when (3^2. Indeed, 
we have 



2 / N 



{2Tizy {\n{TTzr) + Cf + z'{e''"^f 



(4.21) 



Remembering that = z(e '^), we recover a known result aX (3 = 2. The correlation 
function n\j^{r) does not have a finite limit when /? = 2, however the truncated correlation 
function n'^^ {r) = ^^^(r) — (?t.+)^, does have a finite limit at the collapse point (3 = 2 

nf+{r) (277^2)2 (\n{nzr) + C)^ , 



P = 2 



(4.22) 



Furthermore, the density correlation functions of the two-component plasma at /3 = 2 have 



been computed exactly in Refs. 



nf_^ir) = -i2nz'YiKoi27rzr)) 



(4.23) 



where Kq is the modified Bessel function of order 0. Using the well-known expansion of 
Kq for small-argument Q| we verify that the small-r asymptotics ()4.22j) are in complete 
agreement with the exact expression ()4.23|) . 

V. THE OPPOSITE CHARGE CASE: Qi > AND Q2 < 

If Qi and Q2 have different signs, say Qi > and Q2 < 0, then the dominant behavior of 
Eq^Q^ir) is always given by ()3.16|) . since we always have (3\Qi + Q2\ < 2, because in order to 
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satisfy the stability condition < (3Qi < 2 and —2 < f3Q2 < 0. The condition for the term 
n = —1 (or n = +1) to become dominant, exposed in Sec. II VI is never satisfied. Therefore, 
in this section, we will be interested in the next subdominant contribution to EQ^Q^{r) in 
some particular cases. 



A. Effective mean-force potential between two opposite external charges 

In this section we are interested in the case Qi = Q > and Q2 = —Q- Notice that in this 
case in the general small-r expansion (j3.15|) . the term n > and —n are of the same order: 
j"-/3Q^+2n(i-/3/4)+n2/3/2_ Furthermore besides the factor r~^'^^ which is present for all n, the 
rest of the power of r is independent of Q. Thus to study the first subdominant contribution 
we need to take into account both terms n = 1 and n = — 1. Actually, a closer look to 
the expansion (|3.15p shows that in the limit Qi — > —Q2, the coefficients of the expansion 
diverge. But the coefficient for n have a first order pole with opposite residue than the one 
for — n, thus adding both contributions gives a finite result. Expanding the coefficients for 
n = — 1 and n = 1 in ()3.15p when Qi + Q2 ^ and adding both contributions gives 



4 



-/^(l) + 41n 



-fib'' 



2 

4-/3 



4 + 4C.,(-f)+.(f) + ,(l-f) + ,(l.f; 



z' 



+ j^2 r-'''^'-' JW, -f3Q, f3) + 0{r-^^'^') (5.1) 



We have also written explicitly the second subdominant contribution of order r~^'^+^~^ which 
actually does not come from the n = —2 and n = 2 terms but from the first subdominant 
contribution to the n = term: the last term in equation ()3.15|) . An explicit expression for 
J{PQ, —PQ,P) is given in the appendix iBl formulas ()B2IIB6|) . 

The next subdominant contribution 0{r~^'^^~^^), not written here explicitly comes from 
a contribution of the descendant field ((90)^(90)^) 
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B. Counter ion density profile around an single external charge Q > 



The density of negative particles (counterions) of the plasma around a single guest charge 
Q > is given by 

--(-) = ^ ^ (e-^,,) • (5-2) 
The leading term in the small-r expansion is of order r"^*^, as explained earlier. Using the 
general expansion (j3.15|) we find that next subleading contributions comes from the n = 1 
term if Q < 1 and from the n = —1 term if Q > 1. Thus, 

(5.3) 

Notice that if Q = 1 both subleading contributions presented above are of the same order, 
in r^~^, and they should be added as explained in the preceding section. The case Q = 1 
will be studied in the next section. 

The limit (3Q — > 2 is quite unnatural here. Notice that when j3Q — > 2, the denominator 
of the leading term in r~^'^ diverges. Then this term will disappear when [3Q = 2. On 
the other hand the next subleading term, from ()5.3|) . case Q > 1, has a finite limit when 
PQ 2. Indeed, in the numerator ji{—PQ, P, P) has a pole of order 1 at /? = 2, but so does 
^gi6Q0^ in the denominator. Thus, we find 

".M~/|-^^^^-. « = 2. (5.4) 

We call this limit "unnatural" because the power of r changes in a non-continuous manner 
from —PQ to P — 2, contrary to the "natural" limits found in Sec. IIVI Also the mechanism 
for finding a finite result is not the same as in Sec. IIVI where we added terms of the same 
order and their poles canceled each other. 

Let us recall the fact that at PQ = 2, strictly speaking we need to add a short-distance 
regularization to the potential created by the external charge, for example by considering 
that the external charge is not a point but a small impenetrable disk. With this regular- 
ization the denominator in the leading term of ()5.3|) will not vanish when PQ = 2 and 
the density n^{r) will still behave as r~^^, but with a prefactor which is not universal, in 
the sense that it depends on the regularization procedure (it depends on the radius of the 
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external disk particle). This is opposed to the situation for n+(r) that has the universal 
limit ()4.1U|) when (3Q = 2, independent of the regularization procedure (provided the radius 
of the external particle is small enough). 



(2) 

C. Density correlation function n_^_{r) 

Putting Q = 1 in equation (jS.lj) allow us to find the small-r expansion of the density 
correlation function n\_{r) between a positive charge of the plasma and a negative one, 



r->0 



^^^(l) + 41n 



4 



nz \ 4-/3 



7(62 



4 + 4C + ^{^) + ^{-^) + ^(1 - ^) + ^(1 + ^) 



+ zV-^J(/3,-/3,/?) + 0(0 I (5.5) 



An explicit expression for the coefficient J(/5, — /?, P) is given in the appendix |Bl see equa- 
tions (IB7«B9jl . 

As a test, it is interesting to study the limit (3^2 and compare it with the exact 



expression from Refs. |23 . 



24| 



njl^(r) = {2nz'){K^{2nzr)Y , /3 = 2 



(5.6) 



where Ki{x) is the modified Bessel of order 1. 

Notice that when /3 — > 2 the terms in r^"^ and r'^~'^^ from equation ()5.5|1 become of the 
same order. Using the expansions ()A9|) for (e''"^) and ()A12|) for /^(l) from appendix |X1 and 
the expansion ()B2njl from appendix iBl for J(/3, — /3,/5) when /5 ^ 2, we find 



nf_{r) = z\-^ + {2'Kz^f- 



+ C + ln(7rz) + 2(ln(7rz) + C)H + ©(r^) + 0(/3 - 2) . (5.7) 



Using the expansion ()A10|) of (e*^'^)^ we can write 



2\2 



z^r-^ + (27r22) 



C + ln(7rzr) 
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We recover the well-known fact that although n\__{r) does not have a finite limit when 
/3 — s> 2, the truncated correlation function rS^]^ {r) = n^^\{r) — n^n^ = n^l(r) — z^(e*'"^)^ 
has a finite limit when /? = 2, 

+ 0(r2), p = 2. (5.9) 

With the known expansion of Ki{x) for small- argument, we verify that the small-r expan- 
sion ()5.9j) is in complete agreement with the exact expression ()5.()|1 . 



(2)T/ N 

nil r 



2; r 



+ {27CZ 



2\2 



-- + C + hi{TTzr) 



VI. CONCLUSION AND PERSPECTIVES 



Using the operator product expansion we have determined the short-distance behavior 
of the effective potential between two charges Qi and Q2 immersed in a two-dimensional 
charge-symmetric two-component plasma. We were also able to determine the short- distance 
behavior of the density profiles around a single external charge Q in the plasma. The colon 
density profile exhibits a change of behavior at (3Q = 2 — (3, from a dependence in for 
(3Q < 2 — /3 to a dependence in r^~^ for jSQ > 2 — j3. On the other hand the leading 
behavior of the counterion density profile is always r~^^ . We have also determined the first 
subleading contribution to the counterion density profile. Specializing to the case (5 = 1, 
we obtained the short-distance behavior of the density correlation functions. 

This work complements those of Refs. 0, ^| in which the large-distance behavior 
of the correlations and related quantities were determined. Although there is not (yet) a 
closed analytical expression for the correlations functions of the two-component plasma, we 
now have both its large- and short-distance asymptotics. 

A possible interesting application of this work is the following. As we mentioned in the 
introduction, the bulk thermodynamic properties of the two-component plasma have been 
determined for (3 < 2 for a system of point particles However, Kalinay and Samaj 2^ 
have devised a method to obtain the thermodynamic quantities of the two-component plasma 
beyond /? = 2 up to /3 < 3 for a system where the charges are small disks of radius a, in the 
limit no"^ — > (n is the density). A basic ingredient of their method is the knowledge of the 
short- distance behavior of the correlation functions in the region 2 < /5 < 4. Their work is 
based on the leading behavior, in r~^, of n^_{r). If the operator product expansion can be 
used beyond /? = 2, a simple generalization of the analysis of Sees. IIV CI and IV CI gives, for 
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2 < /5 < 4, 

=, - '(^''^ 7(1 - f + f ) + 0(r-^ r«-^^) (6.1) 

n^'l(r) = + A"-2^J(/3,-/3,/3) + 0(r2-^,r«-3'^). (6.2) 

With this expansion, using the method of Ref. the thermodynamic properties of the 
two-component plasma, in the low-density limit, can be obtained up to /? < 10/3. 
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APPENDIX A: EXPANSIONS OF {e'^Q^) 



In the text we need the expansion of (e '^'^) around a regular value of Q and also around 
the pole l3Q = 2. If jSQ < 2, the expansion is immediate to obtain 



ibiQ+e)<f>\ ^ /JbQ<t>\ , I , O I , 

e^O^ dQ ^2 dQ^ 



Using the explicit expression ()3.3|) for (e*^'^'^) we have 



+ 0(e^). 



(Al) 



hiQ) + In ■ 



4-/5 7(62)_ 



+ 2 



2pQ , nz 



4-/3 7(62) 



2p , 7r;z 



+ ^6'(Q) + T^ln 



where 



dhiQ) ^ , I ^ 
dQ 



^ dt 

T 



4 - (3 7(&') 



tsinh(462Qt) 



+ 0(e' 



sinhtcosh[(l — b'^)t\ sinh(62t) 



and 



62^2 cosh(462Qt) 



sinhtcosh[(l — b'^)t\ sinh(62t) 



(A2) 



(A3) 



(A4) 
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If 00 — > 2, the expansion of (e '^'^) can be done following the same steps of the appendix 
of Ref. lOf. Let /SQ = 2 — /Se. First, we isolate the pole by writing Ib{Q) = ^i(e) + ^2(e) + 
J3(e), with 

2 



Xi(e) 



I,(e) 



^ dt 

't 

°° dt 

1 ~t 



2sinh(62t) sinhtsinh[(l - b^)t] 



2 1^5 -el .V- 



2sinh(62t) sinhtsinh[(l - b^)t] 



21 ^-el .V- 



(A5a) 
(A5b) 



dt = r(0, = - ln(46'e) -C + A¥e - 4(6'e)" + O(e^) (A5c) 



where r(a, x) = e^*t"^^ dt is the incomplete Gamma function. The integral X3 contains 
the pole at e = 0, while the integrals Ii and I2 are regular at e = 0. The expansion 
of Ti and X2 can simply be obtained with their Taylor series around e = 0: Xi^2(e) = 
Ji,2(0) + eJi,2(0) + (eV2)Xr,2(0) + O(e3). 

Putting together all the contributions from Xi, X2, X3 and from the expansion of 
(7r;z/7(62))Q'? in (jHSj), we finally find 



1 ^ i?i(/5/4) 



4 , irz 
m ■ 



2-/3g 



^2(/5/4) + -^ln^^ + lf/?i(/?/4)--^ln. 



4-/3 7(/?/4) 2 



4-/3"^7(/?/4)) 



+0((2-/3g)^ 



(A6) 



with 



Ri{b') 



+ t6M 2 



cosht 



26M f 



sinh(62t) cosh[(l - 

4 I - e 



,99/ cosh t coth t + sinh t . 9, 

52^2 / ^ \ ^-2t 



(A7) 
(A8) 



^sinh(62t) cosh[(l - 

Actually, in the main text, we only use the expansion up to order 0(1), the 0{(3Q — 2) term 
is given here only for informational purposes. 

It is also useful to know the expansion of (e*'"^) when (3^2. From the explicit for- 
mula ()2.7p we find 

2txz 



/3-2 2 - /3 



'27Tz[C + \n{TTz)]-7Tz{f3 -2)[C + \n{7Tz)][l + C + \n{TTz)] + 0{{(3 -2y) (A9) 
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and 

(AlOl 

In Sec. IV (^I we need the expansion of /^(l) when (3^2 which is a pole of first order. To 
isolate the pole contribution, it is convenient to write /^(l) as 



4e- 



-2t 



+ 



sinh(/3t) 



t ' sinhtsinh(/5t/4) cosh[(l - f )t] 



- 4e-(2-/3)i 



dt. (All) 



The remaining integral is now convergent when /5 ^ 2, and can be easily expanded around 
;3 = 2. We find 

4 



+ AC - A{(3 - 2) ln2 + 0((/? - 2)^) . 



(A12) 



APPENDIX B: COMPUTATION OF J(/3Q, -/3g, (3) 



In this appendix we give an explicit expression for the integral 



J{(3Q,-[3Q,(3) = / d'xd'y 





X 


|/3Q 


1 


-y 


|/3Q 


\y 




1 - 


X 




x-y 
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(Bl) 

In the appendix D of Ref. |20[, an integral more general than (IBlj) . with arbitrary powers in 
each factor is computed. Using the general results from Ref. |20i], we find for our particular 
case that 

J(/3Q, -/?Q,/3) = s{(3/2f {(Jr)2 + (J-)2} (B2) 

where 



Jr = -^^(^w/2)Ji+ + ^(/?(Q + i)/2)4) 

J2 = -^P(s(/5(Q-l)/2)J+ + .(/3Q/2)J+) 



where we used the notation s(x) = sinfTra;) and 
^1^ 



r(i-f)r(2-f)r(i-^)2 



r(2 _ /3(^)r(3 - ^ 



X 



32I 2'2' 2' 2' 2' 



r(i-f)r(2-f)r(i + ^)^ 
r(2 + ^)r(3 + ^; 



X 



3i^2(2--,— ,1 + — ;2 + ^ ,3 + ;1 



(B3) 
(B4) 



(B5) 



(B6) 
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where 3^2(01, 02, 03; 61, 62; -s^) = YlT=o ^"^{bi)k{b^)kk^'' ^ generalized hypergeometric func- 
tion, and (a)jfc = Via + A;)/r(a) is the Pochhammer symbol. 

For the particular case (5 = 1, the above expressions are simphfied: 



J{P,-(3,(3) = [s{(3/2)J+ + s{P)4Y+[s{P/2)J+]'- 



(B7) 



with 



+ 



Jt 



r(i-f)3r(2-f) 
r(2-/?)r(3-/3) 



F2(l-^,2-^,-^;2-/3,3-/3;l 



r(i-f)r(2-f)r(i + f)2 ^ 



3F2(2--,l + -,-;2,3;l) . 



(B8) 
(B9) 



To study the hmit /3 2, we need the expansion of J(/3, —(3,13) around (3 = 2. In the 
hypergeometric function defining J2 i (3 — 2 is & regular point. We can obtain its expansion 
with the Taylor series 



3F2 ( 2 - ^, 1 + ^, ^; 2, 3; 1 ) 3^2(2, 1, 0; 2, 3; 1) 



2' 2 



+ (/3-2) 



a3F2(2-f,l + f,f;2,3;l) 



d(3 



+ 0{{l3-2f). (BIO) 



p=2 



Deriving under the summation we find 

/3 1 , /3 P. 



a3F2(2-f,l + f,f;2,3;l) 



dl3 



-l+C + t({2 + k) _Ti^_^ 



,.2 S (^ + l)(fc + 2) 6 



(Bll) 



Then 



3F2(2-|l + H2,3;l]^=^2 + (/3-2) - l] + 0{{P - 2f) . 



(B12) 



s{P/2)Jt ^ + ^ (6 + (/3 - 2) + - 2)2) 



With this result we have the following expansions, needed for the expansion of J(/9, —(3, (3), 

(B13) 
(B14) 



s{P)4 =^ -27r - ^ (6 + n') {P - 2) + 0{{(3 - 2f) 



For the expansion of the hypergeometric function in , it is convenient to write it as 



3F2(l-^,2-^,-^;2-/3,3-/3;l'l=l-i^4! + (l-fW) 



2' 2 



2(3-/3) 



(B15) 
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with 



Sill) 



E 

k=2 



(/j/2)(2-f),_,(2-f),_i(2-f), 
2(3-/3)fc_i(3-/3)fcfc! 



(B16) 



This last term can be expanded to the first order using S{P) = S'(2) + (/9— 2)S"(2)+0(/9— 2), 
with S'{2) = - i+c-4>{k-r)+m+^ik+i) ^ _(i + vrVl2)/2. Then, 



S{P) 



4k{k-l) 
/3^2~2 2~ 



TC 



l + ^]+0{{f3-2r 



(B17) 



and finally 

2'' 



3F2(l-^,2-^,-^;2-A3-/?;l 



^'-^-^+^-^(---l)+0iiP-2f). 
/3^2 2 4 8 V 6 ! ' 

(B18) 

We expanded up to order (/? — 2)^, since J^*" has a pole of order 2 at /3 = 2, as opposed to 
J 2 which as a pole of order 1. With this expansion ()B18j) . we find 

2tx 



s{(3/2)Ji 



TT 



+ TT + -(3 + 7r2)(/? - 2) + 0((/? - 2)^) 

/3->2 p — 2 D 



Putting together equations ()B13|) . ()B14|) . and ()B19|) we finally find 

1 1 



J(A -/3,/3) = (2n)' 



(/5-2)2 /3-2 



+ 0(/?-2). 



(B19) 



(B20) 
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